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Summary. Various multi-agent models of wealth distributions defined by micro- 
scopic laws regulating the trades, with or without a saving criterion, are reviewed. 
We discuss and clarify the equilibrium properties of the model with constant global 
saving propensity, resulting in Gamma distributions, and their equivalence to the 
Maxwell-Boltzmann kinetic energy distribution for a system of molecules in an ef- 
fective number of dimensions Da, related to the saving propensity A [M. Patriarca, 
A. Chakraborti, and K. Kaski, Phys. Rev. E 70 (2004) 016104]. We use these re- 
sults to analyze the model in which the individual saving propensities of the agents 
are quenched random variables, and the tail of the equilibrium wealth distribution 
exhibits a Pareto law f(x) oc a; - "" 1 with an exponent a = 1 [A. Chatterjee, B. K. 
Chakrabarti, and S. S. Manna, Physica Scripta T106 (2003) 367]. Here, we show that 
the observed Pareto power law can be explained as arising from the overlap of the 
Maxwell-Boltzmann distributions associated to the various agents, which reach an 
equilibrium state characterized by their individual Gamma distributions. We also 
consider the influence of different types of saving propensity distributions on the 
equilibrium state. 



1 Introduction 

'A rich man is nothing but a poor man with money' — W. C. Fields. 

If money makes the difference in this world, then it is perhaps wise to 
dwell on what money, wealth and income are, to study models for predicting 
the respective distributions, how they are divided among the population of a 
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given country and among different countries. The most common definition of 
money suggests that money is the "Commodity accepted by general consent 
as medium of economics exchange" [1]. In fact money circulates from one eco- 
nomic agent (which can be an individual, firm, country, etc.) to another, thus 
facilitating trade. It is "something which all other goods or services are traded 
for" (for details sec Rcf. [2]) . Throughout history various commodities have 
been used as money, termed usually as "commodity money" which include rare 
seashells or beads, and cattle (such as cows in India). Since the 17th century 
the most common forms have been metal coins, paper notes, and book-keeping 
entries. However, this is not the only important point about money. It is worth 
recalling the four functions of money according to standard economic theory: 

(i) to serve as a medium of exchange universally accepted in trade for goods 
and services 

(ii) to act as a measure of value, making possible the determination of the 
prices and the calculation of costs, or profit and loss 

(iii) to serve as a standard of deferred payments, i.e., a tool for the payment 
of debt or the unit in which loans are made and future transactions are 
fixed 

(iv) to serve as a means of storing wealth not immediately required for use. 

A main feature that emerges from these properties and that is relevant from 
the point of view of the present investigation is that money is the medium in 
which prices or values of all commodities as well as costs, profits, and trans- 
actions can be determined or expressed. As for the wealth, it usually refers 
to things that have economic utility (monetary value or value of exchange), 
or material goods or property. It also represents the abundance of objects 
of value (or riches) and the state of having accumulated these objects. For 
our purpose, it is important to bear in mind that wealth can be measured in 
terms of money. Finally, income is defined as "The amount of money or its 
equivalent received during a period of time in exchange for labor or services, 
from the sale of goods or property, or as profit from financial investments" 
[3] . Therefore, it is also a quantity which can be measured in terms of money 
(per unit time). Thus, money has a two- fold fundamental role, as (i) an ex- 
change medium in economic transactions, and (ii) a unit of measure which 
allows one to quantify (movements of) any type of goods which would other- 
wise be ambiguous to estimate. The similarity with e.g., thermal energy (and 
thermal energy units) in physics is to be noticed. In fact, the description of 
the mutual transformations of apparently different forms of energy, such as 
heat, potential and kinetic energy, is made possible by the recognition of their 
equivalence and the corresponding use of a same unit. And it so happens that 
this same unit is also the traditional unit used for one of the forms of energy. 
For example, one could measure energy in all its forms, as actually done in 
some fields of physics, in degree Kelvin. Without the possibility of expressing 
different goods in terms of the same unit of measure, there simply would not 
be any quantitative approach to economy models, just as there would be no 
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quantitative description of the transformation of the heat in motion and vice 
versa, without a common energy unit. 

2 Multi-agent models for the distribution of wealth 

In recent years several works have considered multi-agent models of a closed 
economy [4, 5, 6, 7, 8, 9, 10, 11, 12]. Despite their simplicity, these models 
predict a realistic shape of the wealth distribution, both in the low income 
part, usually described by a Boltzmann (exponential) distribution, as well 
in the tail, where a power law was observed a century ago by the Italian 
social economist Parcto [13]: the wealth of individuals in a stable economy 
follows the distribution, F(x) oc where F(x) is the upper cumulative 

distribution function, that is the number of people having wealth greater than 
or equal to x, and a is an exponent (known as the Pareto exponent) estimated 
to be between 1 and 2. In such models, N agents exchange a quantity x, 
that has sometimes been defined as wealth and other times as money. As 
explained in the introduction, here money must be interpreted all the goods 
that constitute the agents' wealth expressed in the same currency. To avoid 
confusion, in the following we will use only the term wealth. The states of 
agents are characterized by the wealths {x n }, n = 1,2, . . . ,N. The evolution 
of the system is then carried out according to a prescription, which defines 
a "trading rule" between agents. The evolution can be interpreted both as 
an actual time evolution or a Monte Carlo optimization procedure, aimed at 
finding the equilibrium distribution. At every time step two agents i and j 
are extracted at random and an amount of wealth Ax is exchanged between 
them, 

x'j = Xj + Ax. (1) 

It can be noticed that in this way the quantity x is conserved during the single 
transactions, x\ + x'j = Xi + Xj. Here x\ and x'j are the agent wealths after the 
transaction has taken place. Several rules have been studied for the model. 

2.1 Basic model without saving: Boltzmann distribution 

In the first version of the model, so far unnoticed in later literature, the money 
difference Ax is assumed to have a constant value [4, 5, 6], 

Ax = Ax ■ (2) 

This rule, together with the constraint that transactions can take place only if 
x[ > and x'j > 0, provides a Boltzmann distribution; see the curve for A = 
in Fig. 1 . An equilibrium distribution with exponential tail is also obtained if 



4 Patriarca, Chakraborti, Kaski, Germano 

Ax is a random fraction e of the wealth of one of the two agents, Ax — exi or 
Ax = exj . A trading rule based on the random redistribution of the sum of the 
wealths of the two agents has been introduced by Dragulescu and Yakovenko 
[7], 

x i "f~ x j ) > 

x'j = e(xi +xj), (3) 

where e is a random number uniformly distributed between and 1 and e is 
the complementary fraction, i.e. e + e = 1. Equations (3) are easily shown to 
correspond to the trading rule (1), with 

Ax = exi — exj . (4) 

In the following we will concentrate on the latter version of the model and 
its generalizations, though both the versions of the basic model defined by 
Eqs. (2) or (4) lead to the Boltzmann distribution, 

/(X)= (^ CXP (~M) 

where the effective temperature of the system is just the tlth (x) 

[4, 5, 6, 7]. The result (5) is found to be robust, in that it is largely inde- 
pendent of various factors. Namely, it is obtained for the various forms of Ax 
mentioned above, for pairwise as well as multi-agent interactions, for arbi- 
trary initial conditions [8], and finally, for random or consecutive extraction 
of the interacting agents. The Boltzmann distribution thus obtained has been 
sometimes referred to as an "unfair distribution" , in that it is characterized 
by a majority of poor agents and very few rich agents, as signaled in particu- 
lar by a zero mode and by the exponential tail. The form of distribution (5) 
will be referred to as the Boltzmann distribution and is also known as Gibbs 
distribution. 



(5) 



2.2 Minimum investment model without saving 

Despite the Boltzmann distribution is robust respect to the variation of several 
parameters, the way it depends on the details of the trading rule is subtle. 
For instance, in the model studied in Ref. [9] , the equilibrium distribution can 
have a very different shape. In that model it is assumed that both economic 
agents i and j invest the same amount x m i n , which is taken as the minimum 
wealth of the two agents, x m i n = min{a:j, Xj}. The wealths after the trade are 
x\ — Xi + Ax and x'j — Xj — Ax, where Ax = (2e — l)x min . 

We note that once an agent has lost all his wealth, he is unable to trade 
because x m i n has become zero. Thus, a trader is effectively driven out of the 
market once he loses all his wealth. In this way, after a sufficient number of 
transactions only one trader survives in the market with the entire amount of 
wealth, whereas the rest of the traders have zero wealth. 
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Fig. 1. Probability density for wealth x. The curve for A = is the Boltzmann 
function f(x) = (x)' 1 exp(~x / (x)) for the basic model of Sec. 2.1. The other curves 
correspond to a global saving propensity A > 0, see Sec. 2.3. 



2.3 Model with constant global saving propensity: 
Gamma distribution 

A step toward generalizing the basic model and making it more realistic is 
the introduction of a saving criterion regulating the trading dynamics. This 
can be achieved defining a saving propensity < A < 1, that represents the 
fraction of wealth saved — and not reshuffled — during a transaction. The 
dynamics of the model is as follows [8, 9]: 

x' t = Xxi + e(l - X)(x l + Xj) , 

x'j = Xxj + e(l — X)(xi + xj) , (6) 

with e = 1 — e, corresponding to a Ax in Eq. (1) given by 

Ax = (1 — X)[exi — exj] . (7) 

This model leads to a qualitatively different equilibrium distribution. In par- 
ticular, it has a mode x m > and a zero limit for small x, i.e. f(x—* 0)— > 0, 
see Fig. 1. The functional form of such a distribution has been conjectured to 
be a Gamma distribution on the base of an analogy with the kinetic theory 
of gases, which is consistent with the excellent fitting provided to numerical 
data [14, 15]. Its form can be conveniently written by denning the effective 
dimension D\ as [15] 

Dx_ } . 3A _1 + 2A 

According to the equipartition theorem, one can introduce a corresponding 
temperature defined by the relation (x) — D\T\/2, i.e. 



(9) 
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Then the distribution for the reduced variable £ = x/T\ reads 

fi0 = nkm ^ /2_1 cxp( ~ e) = 1Dj2[i) ' (10) 

i.e. a Gamma distribution of order D\/2. For integer or half-integer values of 
n = -Da/2, this function is identical to the equilibrium Maxwcll-Boltzmann 
distribution of the kinetic energy for a system of molecules in thermal equi- 
librium at temperature T\ in a Z?A-dimensional space (see Appendix A for 
a detailed derivation). For D\ = 2, the Gamma distribution reduces to the 
Boltzmann distribution. This extension of the equivalence between kinetic the- 
ory and closed economy models to values < A < 1 is summarized in Table 
1. This equivalence between a multi-agent system with a saving propensity 



Table 1. Analogy between kinetic and multi-agent model 



Kinetic model 



Economic model 



variable 
units 

interaction 
dimension 
temperature 
reduced variable 



K (kinetic energy) x (wealth) 
N particles N agents 

collisions 
integer D 
k B T = 2 (K) /D 
i = K/k B T 



trades 

real number D\ [see Eq. (8)] 
T\ = 2 (x) /D x 



equilibrium distribution /(£) = 7 D / 2 (£) /(£) = 7d a /2(£) 



< A < 1 and an iV-particle system in a space with effective dimension 
D\ > 2 was originally suggested by simple considerations about the kinetics 
of a collision between two molecules. In fact, for kincmatical reasons during 
such an event only a fraction of the total kinetic energy can be exchanged. 
Such a fraction is of the order of 1 — A w to be compared with the 

expression 1 — A = 3/(D/2 + 2) derived from Eq. (8) [15]. While A varies 
between and 1, the parameter D\ monotonously increases from 2 to oo, and 
the effective temperature T\ correspondingly decreases from (x) to zero; see 
Fig. 2. It is to be noticed that according to the equipartition theorem only 
in D\ = 2 effective dimensions (A = 0) the temperature coincides with the 
average value (a;), T\ = 2 (x) /2 = (x), as originally found in the basic model 
[4, 5, 6, 7]. In its general meaning, temperature represents rather an estimate 
of the fluctuation of the quantity x around its average value. The equiparti- 
tion theorem always gives a temperature smaller than the average value (x) 
for a number of dimensions larger than two. In the present case, Eqs. (8) or 
(9) show that this happens for any A > 0. 

The dependence of the fluctuations of the quantity x on the saving propen- 
sity A was studied in Ref. [8]. In particular, the decrease in the amplitude of 
the fluctuations with increasing A is shown in Fig. 3. 
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Fig. 3. Reproduced from Ref. [8] (only here m = x while T = 1 is a constant). The 
continuous and the dotted curves are the wealths of two agents with A = 0.9 and 
A = 0.5 respectively: notice the larger fluctuations in correspondence of the smaller 
A. The inset shows that Am = Ax = y ((a; — (x)) 2 ) decreases with A. 



The fact that in general the market temperature T\ decreases with A means 
smaller fluctuations of x during trades, consistently with the saving criterion, 
i.e. with a A > 0. One can notice that in fact T\ = (1 — A) (x) /(l + 2A) ss 
(1 — A) (a;) is of the order of the average amount of wealth exchanged during 
a single interaction between agents, see Eqs. (6). 
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2.4 Model with individual saving propensities: Pareto tail 

In order to take into account the natural diversity between various agents, 
a model with individual propensities {Xi} as quenched random variables was 
studied in Refs. [10, 11]. The dynamics of this model is the following: 



where, as above, e = 1 — e. This corresponds to a Ax in Eq. (1) given by 



Besides the use of this trading rule, a further prescription is given in the 
model, namely an average over the initial random assignment of the individual 
saving propensities: With a given configuration {A^}, the system is evolved 
until equilibrium is reached, then a new set of random saving propensities {A^} 
is extracted and reassigned to all agents, and the whole procedure is repeated 
many times. As a result of the average over the equilibrium distributions 
corresponding to the various {A^} configurations, one obtains a distribution 
with a power law tail, f(x) oc x~ a ~ 1 , where the Pareto exponent has the value 
a = 1. This value of the exponent has been predicted by various theoretical 
approaches to the modeling of multi-agent systems [16, 17, 12]. 

3 Further analysis of the model with individual saving 
propensities 

On one hand, the model with individual saving propensities relaxes toward a 
power law distribution — with the prescription mentioned above to average 
the distribution over many equilibrium states corresponding to different con- 
figurations {Xi}. On the other hand, the models with a global saving propen- 
sity A > and the basic model with A = 0, despite being particular cases 
of the general model with individual saving propensities, relax toward very 
different distributions, namely a Gamma and a Boltzmann distribution, re- 
spectively. In this section we show that this difference can be reconciled by 
illustrating how the observed power law is due to the superposition of differ- 
ent distributions with exponential tails corresponding to subsystems of agents 
with the same value of A. 

3.1 The x-\ correlation 

A key point which explains many of the features of the model and of the 
corresponding equilibrium state is a well-defined correlation between average 
wealth and saving propensity, which has been unnoticed so far in the literature 



x\ = XtXi + e[(l - Xi)Xi + (1 - Xj)xj] , 
x'j = XjXj + e[(l - Xi)xi + (1 - Xj)xj] , 



(11) 



Ax = e(l — Xi)xi — e(l — Xj)xj . 



(12) 
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Fig. 4. Equilibrium state in the x-X plane after t — 10 trades for a system of 
TV = 1000 agents. Left: Circles (o) represent agents, crosses (x) represent the average 
wealth (x(A)}, the continuous line is the function (x(X)) = k/(1 — A), with k — 
0.1128. Right: The product (x(X)) (1 — A) (dotted line) is constant, in agreement 
with Eq. (14). The product Ax(X)(l — A) (continuous line), where Zia;(A) is the 
standard deviation, shows that Ax(X) grows slower than (x(X)). 



[18]. The existence of such a correlation can be related to the origin of the 
power law and its cut-off at high values of x. It also explains the paradox 
according to which a very rich agent may lose all his wealth when interacting 
with poor agents, as a consequence of the stochastic character of the trading 
rule defined by Eq. (11). Figure 4 shows the equilibrium state for a system with 
N = 1000 agents after t = 10 9 trades. Each agent is represented by a circle (o) 
in the wealth-saving propensity Xi-Xi plane. The correlation between wealth x 
and saving propensity A becomes very high at large values of x and A. Namely, 
one observes that the average wealth (x(X)} [crosses (x) in Fig. 4] diverges 
for A —* 1. The average (x(A)) was obtained by computing the probability 
density f(x,X) in the x-X plane (normalized so that J dx dX f(x, A) = 1) and 
averaging for a fixed value of A, 

(x(X)) = J dxxf(x, A) . (13) 

The observed correlation naturally follows from the structure of the trade 
dynamics (11). We remind that initially every agent has the same wealth 
£0 = {x}. During the initial phase of the evolution, when all agents have 
approximately the same wealth (x) , an agent i with a large saving propensity 
Xi can save more — on average — and therefore accumulate more. Afterwards, 
the agent i will continue to enter trades by investing only a small fraction 
(1 — Xi)Xi of his wealth Xi in the trade. Even when interacting with an agent 
j, with a smaller wealth Xj < Xi, agent i may still be successful in the trading, 
since agent j may have also a smaller saving propensity Xj , so that the traded 
fraction of wealth (1 — Xj)xj is comparable with or even larger than (1 — X{)Xi. 
Trading by agent i will very probably be successful on average with all agents 
j with a Xj such that (l-Xj)xj is smaller than (1 — Aj)xj. These considerations 
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Fig. 5. Histogram of the wealths x' = x(t + 1) after a trade versus x — x(t) before 
the trade for all agents and trades, in a system with N = 1000 agents and 10 9 trades. 
The distribution is narrower for large x (rich agents), implying that it is unlikely 
that a rich agent becomes poor within a single trade. 



suggest that agent i will reach equilibrium (and his maximum possible wealth) 
when (1 — XijXi = K « ((1 — X)x). The ratio between the constant k and the 
average ((1 — X)x) = $^,-(1 — Xj)xj/N is actually found to be of the order of 
magnitude of 10. The formula 

(14) 

however, shown as a continuous line in Fig. 4, provides an excellent inter- 
polation of the average wealth (x(A)) (also shown in the figure) computed 
numerically. 

3.2 Variation of a single agent's wealth 

The stability of the asymptotic state is also shown by the histogram in Fig. 5 
of the wealths x' = x(t + 1) after a trade versus x = x(t) before it defined 
in Eqs. (1). The distribution is narrower at larger values of x than at smaller 
ones, implying that the probability that an agent i will undergo a large rel- 
ative variation of his wealth Xi within a single trade is much higher for poor 
agents. The situation at small x (corresponding to agents with smaller saving 
propensities) is instead more similar to that of the trading rule without saving 
(A = 0), Eqs. (6): the distribution is broader, indicating a higher probability 
of a large wealth reshuffling during a trade. 

3.3 Power laws at small x and t scales 

A peculiarity of the model with individual saving propensities is noteworthy. 
On one hand, in the procedure used to obtain a power law in Ref. [10] all agents 
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Fig. 6. The equilibrium configurations corresponding to four different random sav- 
ing propensity sets {Xi}, for a system with N = 1000 agents, differ especially at 
large x where the distribution deviates from a power law. 



are equivalent to each other: they enter the dynamical evolution law on an 
equal footing, their saving propensities are reassigned randomly with the same 
uniform distribution between and 1, and even their initial conditions can be 
set to be all equal to each other, Xi — (x), without loss of generality. Therefore 
the various equilibrium configurations, corresponding to different sets {A^}, 
are expected to be statistically equivalent to each other, in the sense that 
one should be able to obtain the power law distribution by a simple ensemble 
average for any fixed configuration of saving propensities {A^}, if the number 
of agents N is large enough. On the other hand, an averaging procedure over 
several {Xi} configurations is in practice necessary to obtain a power law 
distribution. 

In order to understand this apparent paradox, we checked how the equi- 
librium distributions, corresponding to a given set of saving propensities {A^}, 
look like. One finds that every configuration {A^} produces equilibrium dis- 
tributions very different from each other; see Fig. 6 for some examples. The 
structures observed are very different from power laws, with well resolved 
peaks at large x. Only when an average over different {A^} configurations 
is carried out, one obtains a smooth power law with Pareto exponent a = 1. 
These same figures show, however, that for a given configuration {Ai} a power 
law is actually observed at small values of x. Another related interesting fea- 
ture of simulations employing a single saving propensity configuration {Ai} is 
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Fig. 7. Time evolution of the x distribution of a system with N = 1000 agents: 
t = 2 x 10 6 (top left), 3 x 10 6 (top-right), 5 x 10 6 (bottom- left), and 2 x 10 7 trades 
(bottom-right). The distribution looks as a power law at small times, but develops 
into a structured distribution, maintaining a power law shape only at small x. 



that a power law distribution is found only on a limited time scale, while it 
disappears partly or totally at equilibrium. Thus also in the time dimension 
one surprisingly finds a distribution much more similar to a power law at a 
smaller rather than larger scale. This is shown in the example in Fig. 7, where 
the distributions of a system of 1000 agents at four different times are com- 
pared to each other. These features suggest that the power law is intrinsically 
built into the dynamical laws of the model but that, for some reasons, it fades 
away at large x and t scales. The x-X correlation discussed above in Sec. 3.1 
can provide an explanation of these features, both for those in the x and in 
the time dimension, as discussed below. 

3.4 Origin of the power law 

The peculiar features illustrated above, the necessity of averaging over differ- 
ent configurations {A^} as done in Ref. [10] to obtain a power law distribu- 
tion, as well as the power law distribution itself, are here explained in terms of 
equilibrium states of suitably defined subsystems and the x-X correlation illus- 
trated above. This may seem odd since at first sight the averaging procedure 
of Ref. [10] defines a nonequilibrium process, the system being brought out of 
equilibrium from time to time by the reassignment of the saving propensities. 
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Fig. 8. Semi-log and log plots of partial distributions (continuous curves) and the 
resulting overlap (dotted line). Above: Partial distributions from the 10 intervals of 
width AX — 0.1 of the total A range (0,1). Below: The last partial distribution (with 
a power law tail) above, corresponding to the interval A = (0.9, 1.0), has been further 
resolved into ten partial distributions for the sub-intervals of width AX = 0.01. 



Correspondingly, one may attribute the power law to the underlying dynam- 
ical process, as it is often the case in nonequilibrium models (e.g. models of 
markets on networks [19]). 

However, if one considers the partial distributions of agents with a certain 
value of A, one finds an unexpected result. For numerical reasons we con- 
sider the subsets made up of those agents with saving propensity A within a 
window AX around a given value A. Figure 8 (upper row) shows the partial 
distributions (continuous lines) of the ten subsystems obtained by dividing the 
A range (0,1) into ten slices of width AX = 0.1 and average values 0.05,. . .,0.95 
(curves from left to right respectively). Most of the partial distributions have 
an exponential tail, and only when summed up their overlap (dotted line) 
reproduces a power law. It can be noticed that the last partial distribution, 
corresponding to the interval A = (0.9, 1.0), is not of exponential form, but 
rather presents a power law tail, which overlaps with the total distribution 
at large x. However, its power law form is due only to the low resolution in 
A employed. In fact even this partial distribution can in turn be shown to be 
given by the superposition of exponential tails. By increasing the resolution 
in A, i.e. using a smaller interval AX = 0.01 to further resolve the interval 
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A = (0.9, 1.0) into subintervals with average values A = 0.905, . . . ,0.995, one 
obtains the partial distributions shown in the lower row of Fig. 8. It is to be 
noticed that also in this case the last partial distribution corresponding to 
the interval A = (0.99, 1.00) has a power law tail. The procedure can then in 
principle be reiterated to resolve also this partial distribution by increasing 
the resolution in A. 

These facts also explain the origin of the peaks visible at large x in the 
plots in Figs. 6 and 7. Due to the high wealth-saving propensity correlation 
at large values of x, see Fig. 4, these peaks are due to agents with high A. The 
reason why these agents give rise to resolved peaks instead of contributing to 
extending the power law tail is that the partial distributions (i.e. the average 
values) of single agents get farther and farther from each other for A — > 1, 
while the corresponding widths do not grow enough to ensure the overlap 
of the distributions of neighbor agents in A-space. Eventually, each agent (or 
cluster of agents) at high values of x will be resolved as an isolated peak against 
the background of the total distribution. In greater detail, one finds that the 
average value (x(X)) diverges for A — ► 1 as 1/(1 — A), as shown in Fig. 4. 
This implies that also the distance between two generic consecutive agents 
increases: if agents are labeled from i = 1 to i — N in order of increasing 
A (Ai < . . . < A^v) and the A distribution is uniform, then AX = Aj+i — 
A.; = constant. The distance between the average positions of the partial 
distributions of two consecutive agents is, from Eq. (14), 

5(x(X)) = (x(X + AX)) {x(X)) « A\ « , (15) 

where n is a constant. Thus S(x(X)) diverges even faster than (x(X)). At the 
same time, the width of the partial distribution Ax(X), here estimated as 

Ax(X) = \J (x 2 (X)) — (x(X)) 2 , grows slower than (a;(A)), i.e. for A — > 1 the 
ratio Ax(X)/ (x(X)) — > 0; see Fig. 4 (right). The breaking of the power law 
and the appearance of the isolated peaks takes place at a cutoff x c where the 
distance 5(x(X)} between the peaks corresponding to consecutive agents i and 
i + 1 becomes comparable with the peak width Ax(X). 

Also the origin of the peculiarities in the time evolution of the distribution 
function, mentioned in Sec. 3.3, can now be explained easily. In order to reach 
the asymptotic equilibrium state, agents can rely on an income flux which 
is on average proportional to Xi(l — Ai). At the beginning, when all agents 
have the same wealth Xi — xq, this quantity is smaller for agents with a 
larger Aj; and with this smaller flux agents with large Ai have to reach their 
higher asymptotic value ocl/(l — Aj). Asa consequence, the relaxation time 
of an agent is larger for larger A, a result already found in the numerical 
simulations of the multi-agent model with fixed global saving propensity (see 
Fig. 2 in Ref. [8]). Correspondingly, partial distributions of rich agents will 
reach their asymptotic form later (last frame in Fig. 7), while, at intermediate 
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Fig. 9. Wealth distribution of a system of N = 10 6 agents after 10 12 trades: the 
uniform A distribution Eq. (16) produces a smoother distribution with a power law 
shape extending to higher x (left) than for a random A distribution (right). 



times, their distributions will be spread at smaller values of x, contributing 
to smoothen the total distribution (first frame in Fig. 7). 

It is also possible to explain why the averaging procedure of Ref. [10] 
is successful in producing a power law distribution. Averaging over different 
configurations {A^} is equivalent to simulate a very dense distribution in A — 
which has large relaxation time and number of agents — with an affordable 
number of agents and computer time. However, the procedure is not needed in 
principle, since the power law can be obtained also when a single configuration 
with a proper density in A-space is used. 

3.5 Checking different A distributions 

A practical way to avoid the appearance of the peaks at large x and obtain a 
distribution closer to a power law is to increase the density of agents, especially 
at values of A close to 1. In a random extraction of {Xi}, it is natural that 
consecutive values of Aj will not be equally spaced. Even small differences 
will be amplified at large x and will result in the appearance of peaks. A 
deterministic assignment of the A, e.g. a uniform distribution achieved through 
the following assignment, 

A, = A i = 0,JV-l, (16) 

is a uniform distribution of A in the interval [0,1) and will generate a smoother 
distribution of x. The comparison of the results for this distribution with those 
for a random distribution of A is done in Fig. 9 (notice the high value of N). 
In the uniform case not only the power law extends to higher values of x but 
also that the distribution of peaks at large x is globally smoother, in the sense 
that on average the single peaks follow a power law better. 
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4 Conclusions 



We have reviewed some multi-agent models for the distribution of wealth, in 
which wealth is exchanged at random in the presence of saving quantified by 
the saving propensity A. We have shown how a distribution of A generates a 
power law distribution of wealth through the superposition of Gamma distri- 
butions corresponding to particular subsets of agents. The physical picture for 
the model with individual saving propensities is thus more similar to that of 
the model with a constant global saving propensity than it may seem at first 
sight. In fact any subset of agents with the same value of the saving propensity 
A equilibrates in a way similar to agents in the model with a global saving 
propensity i.e. leading to a wealth distribution with an exponential tail. Cor- 
respondingly we have shown that both the noise in the power law tail and 
the cutoff in the power law depend on the coarseness of the A distribution. 
This extends the analogy between economic and gas-like systems beyond the 
case of a global A > 0, characterized by a Maxwcll-Boltzmann distribution, to 
uniform continuous distributions in A that span the whole interval A £ [0, 1). 
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A Maxwell-Boltzmann distribution in D dimensions 

Here we show that for integer or half-integer values of the parameter n the 
Gamma distribution 



where -T(n) is the Gamma function, represents the distribution of the rescaled 
kinetic energy £ = K/T for a classical mechanical system in D = 2n dimen- 
sions. In this section, T represents the absolute temperature of the system 
multiplied by the Boltzmann constant k^. 

We start from a system Hamiltonian of the form 



7n(0 = -T(»r 1 f l - 1 exp(-0, 



(17) 




(18) 
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where P = {pi, . . . , p^} and Q = {q l7 . . . , q N } are the momentum and po- 
sition vectors of the N particles, while V^(Q) is the potential energy con- 
tribution to the total energy. For systems of this type, in which the to- 
tal energy factorizes as a sum of kinetic and potential contributions, the 
normalized probability distribution in momentum space is simply /(P) = 
Yl^irniiT)- / 2 cxp(-p 4 2 /2m l T). Thus, since the kinetic energy distribution 
factorizes as a sum of single particle contributions, the probability density fac- 
torizes as a product of single particle densities, each one of the form 

(19) 



/(P)= (W)^ CXP ('2^t) ' 



where p = (pi, . . . ,pd) is the momentum of a generic particle. It is convenient 
to introduce the momentum modulus p of a particle in D dimensions, 



P 2 -P 2 = f> 2 fe , (2°) 
fe=i 

where the pk's are the Cartesian components, since the distribution (19) de- 
pends only on p = . One can then integrate the distribution over the 
D — 1 angular variables to obtain the momentum modulus distribution func- 
tion, with the help of the formula for the surface of a hypcrsphere of radius p 
in D dimensions, 

"■M-mm*^ 1 - (21) 

One obtains 

,D— 1 ( P' 



HP) = SMtip) = r{Dm2mTrP P-^ (-J-f) ■ P*> 
The corresponding distribution for the kinetic energy K = p 2 /2m is therefore 



f(K) 



d P a \ 



1 (K\ D ' 2 - 1 ( K 



p=V2mK 



r(D/2)T\TJ *V T 



exp [-■=). (23) 



Comparison with the Gamma distribution, Eq. (17), shows that the Maxwell- 
Boltzmann kinetic energy distribution in D dimensions can be expressed as 

f{K)=T- 1 lD/2 {K/T). (24) 

The distribution for the rescaled kinetic energy, 

Z = K/T, (25) 

is just the Gamma distribution of order D/2, 

/(0 



"4 /(*)] ^ - J^J2) ^ exp(-0 = 7^(0 • (26) 
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